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1 . SUMMARY 


A computer program to calculate the radiation properties 
of the reflector antennas is presented. In [ts present form 
can be used for paraboloidal, spherical, or 

tne program 

ellipsoidal reflector surfaces* However, it can be easily 
modified to handle any surface that can be expressed analytically. 
The program is general enough to allow any arbitrary location 
and pointing angle for the feed antenna. The effect of blockage 
due to the feed horn is also included in the computations. 

The computer program is based upon the technique of tracing 
the rays from the feed antenna to the reflector to an aperture 
plane. The far field radiation properties are then calculated 
by performing a double integration over the field points in 
the aperture plane. To facilitate the computation of double 
integral , hij^ver, the field points are first aligned along 
the equispaced straight lines in the aperture plane. The 
computation time is relatively insensitive to the absolute size 
of the aperture and even though no limits on the largest reflector 
size have been determined, the program has been used for 
reflector diameters of 1000 ^. 


2. INTRODUCTION 


The purpose of this report is to present a detailed 
description of a computer program called REFLCTR which has 
been written to calculate the radiation properties of reflec- 
tor antennas whose surfaces can be analytically expressed in 
Cartesian coordinates. The reflector feed can have any 
arbitrary location (does not have to be at focus) and arbitrary 

t 

pointing angle. In its present fora, the program is set up 

2 2 

to handle either of a paraboloidal (y +z *4f (f+x) where f is 

2 2 2 2 

the focal length ) , spherical (x +y +z =r where r is the radius) , 

2 2 2 2 2 

or ellipsoidal ( (x /a )+(y +z )/b =1 where a and b are semi 
major ^and semi minor axes) reflector surface whereas the 
reflector can be defined by the curved surface of a frustum 
of any one of the above mentioned shapes. 

The theory associated with the computer program is briefly 
presented in Section 3 and then the computer program is presented 
in Section 4. An example is presented in Section 5 for a 
spherical reflector to demonstrate the methodology of using 
the computer program. A detailed description of some of the 
subroutines is given in Sections 6-10. Section 11 shows how 
the program can be very easily modified to handle any other 
reflector surface that can be expressed analytically. 


3. THEORY 


The underlying theory is described in detail in 
Reference 1. For the purpose of understanding the computer 
program, however, a brief description of the theory is 
presented . 

Consider a reflector and its feed as shown in Figure 1. 
The reflector surface is known in the (x,y,z) coordinate 
system, 0 being the origin of this reflector coordinate 
system. The feed antenna which is allowed to have any 
arbitrary location and direction is directed along the 
negative x'-axis of its own coordinate system whose origin 
is O'. The radiation pattern of the feed antenna is assumed 
to be known in the primed or the feed coordinate system. 

And therefore, to be able to write the equation of a ray 
emanating from the feed antenna in th-*» reflector coordinate 
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Figure 1 — Feed and Reflector Coordinate Systems 


system, a transformation from the feed to the reflector 

coordinates is needed. In general, these two coordinate 

- » 

systems are related to each other by a translation and 
a rotation, each of which will have three orthogonal com- 
ponents. The computer program is written such that the 
translation is defined as the vector O'O expressed in the 
primed coordinate system i.e . , the coordinates of point 0 
in the feed coordinate system. The three components of 
this vector 0*0 are called FEED(l) , FEED (2) and FEED (3) 
in the computer program. 

The rotation needed to relate the primed and the unprimed 
coordinate systems also has three components which are called 
ALPHA, BETA, and GAMMA in the computer program. These three 
angles are defined such that if the primed coordinate system 
is rotated by ALPHA, BETA, and GAMMA, the primed coordinate 
axes become parallel to the unprimed coordinate axes. Speci- 
fically, ALPHA is the rotation of the x*y' -plane about the 
z 1 -axis needed to make the y*-axis parallel to the y-axis, 
and as shown in Figure 2(a) for a simple case, ALPHA allows 
the 'feed to be pointed anywhere along a horizontal line on 
the reflector. BETA is the rotation of the y'z' -plane about 
the x' -axis to align the z'-axis with the z-axis a> i allows 
the feed antenna to have any arbitrary direction of polarization, 
i.e., not necessarily coincident with either the z- or the 
y-axis as shown in Figure 2(b). Finally, GAMMA is the rotation 
of the x'z* -plane about the y'-axis to make the x'-axis 
parallel to the x-axis. <?AMMA allows the feed to be pointed 



Figure 2 — Angles ALPHA, BETA, and GAMMA 
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anywhere along a vertical line on the reflector as shown in 
Figure 2 (c) . The rotation angles are considered positive 
for counter clockwise rotation when looking in the negative 
direction along the axis of rotation. Thus, for example, if 
the feed antenna was directed upward along the offset focal 
axis by an angle 6 q, GAMMA will be -8 q. To describe the 

j 

general orientation of a feed, however, all three of ALPHA, 
BETA, and GAMMA will be needed. 

The feed pattern is needed at enough equally separated 
values of 6' and so that the rays emanating from the feed 
over these ranges of 8' and more than illuminate the 
reflector, i.e., the rays with the upper and the lower limit 
values of 6* and miss the reflector. A typical ray along 
(8!^, $j) direction is shown in Figure 3. Index i is used to 
denote the change in the angle 8 ' which has NT equispaced 
values between 8^ and 8^. Similarly, index j is used to 
denote the change in the angle which has NP equispaced 
values between ♦ £ and making the total number of rays 
considered emanating from the feed equal to the product of NT 
and NP. As mentioned above, the limits on 6' and are 
required to be such that at least the rays with 6' » 8| and 
8 ^ t and the rays with ♦' * and miss the reflector. 

The reason for this is explained later. 

Associated with each ray (i.e., with each (i,j)) are five 
quantities — Eg and E^, the two components of the electric 
field along the ray, the phase of the electric field, and the 


r 



Figure 3 -- Feed Pattern Angles 

angles (9£,$.J) specifying the direction of the ray. These 
five quantities are stored for each ray in a three-dimensional 
array P(5,NT,NP) which is shown pictorially in Figure 4. 

Using the three pointing angles ALPHA, BETA, and GAMMA of 
the feed with respect to the (x,y,z) coordinate axes, the 
Euler Rotation Matrix (dimensioned as A(3,3)) is calculated. 
Use of the matrix A allows one to write a unit vector along 
each of the rays in the (x,y,z) coordinate system, following 
which, the point of intersection of the ray with the reflector 
(x 0 ,y 0 ,z 0 ) is found, and the x-, y-, z-components of the unit 
normal, NHAT(l), NHAT ( 2 ) , and NHAT ( 3 ) are calculated 
at that point (Figure 5). Knowing the vector along the 
incident ray and along the normal at the point of incidence. 
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the x-, y-> z-components of the vector along the reflected 
ray, SR(1), SR(2) , and SR(3) are calculated. Also, with 
the knowledge of the unit normal vector, and Eg and of 
the incident ray, the three cartesian components of the 
electric field in the reflected ray, ER(1) , ER(2) , and 
ERO) are evaluated. 

A plane parallel to the yz-plane is defined as the 
aperture plane. X is the distance of this plane from the 
origin and it is chosen such that the aperture plane lies 
in front of and near the edge of the reflector. The boundary 
of the aperture plane is a contour formed by the intersection 
of the aperture plane and the rays reflected from the edge 
of the reflector. This contour in general can be a combination 
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Figure 4 — P Array and the Feed Pattern Information 


of two ellipses. For the sake of simplicity in calculation, 
though, the boundary of the aperture plane is approximated 
by a single best fit ellipse. The parameters that define this 
ellipse are its half major axis (HFMAEX, along the y-direction) , 
the half minor axis (HFMXEX, along the z-direction) , and the 
coordinates of the center of the ellipse (XC, YC, ZC) . 

Those rays emanating from the feed which will miss the 
reflector are assumed to be reflected as though the reflector 
surface extended past its boundary. The pcy.nt of intersection 
of each of the reflected rays with the aperture plane called 
(Y,Z) is calculated. Using the points which fall immediately 
inside and outside the aperture plane elliptical boundary, 
the edge points are interpolated at the boundary of the ellipse. 



original 
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Ths points outsids ths aperture plane elliptical boundary 
are then discarded. For th« points that fall inside the 
apfcrture plane ellipse and the ones interpolated on the 
ellipse boundary, Y,Z, ER(2),ER(3), and phase are stored 
in the P array. 

The points of intersection of the reflected rays with 
the aperture plane are in general, n&t uniformly distributed 
over the aperture plane. However, in order to facilitate 
the calculation of the far field radiation pattern from the 
aperture plane distribution, these points in the aperture 
plane are aligned along equispaced y*constant lines. Or, 
equivalently, the magnitude and the phase at each of the 
intersection points are arbitrarily assigned to the nearest 
point on a y«constant line. In this way, the double integral 
to be performed cn the aperture plane reduces to an integral 
over straight line integrals along y«const.ant linos. 
spacing between the y-constant lines must be chosen very 
carefully. For example, shown in Figure 6 is a computer 
generated plot of the points of intersection of the reflected 
rays and the aperture plane. Observe that, as stated earlier, 
these are not uniformly distributed. Figure 7 is another 
computer generated plot showing the locations of the same 
points after they have been aligned alonq reuispaccd y=constnnt 
lines. Notice the blanks created in the lines as a result of 
aligning the points. The spacing between v*const.ant lines 
needs to be chosen to minimize these blanks and to place 
them as near to the aperture edge flow field strength reoi on' 


as possible. The best value for lire spacing has been found 
to be approximately equal to the average spacing of the points 
(Figure 6) along the y-direction. Since the integration time 
depends upon the number of points (of intersection of the re- 
flected ravj and the aperture plane) which in turn is determined 
by the angular increment between the rays used in the feed 
pattern, the integration time is relatively insensitive to the 
absolute size of the reflector. 
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Figure 7 -- Points of Intersection of the Reflected Rays 
and the Aperture Plane after Aligning along 
Equispaced y* Constant Lines 




4. COMPUTER PROGRAM 


The computer program has a main program that calls 
several subroutines, some of which call other subroutines. 
The order in which the subroutines are called is illustrated 
in Figure 8. The remainder of this section describes the 
overall working of the program and the subroutines. The 
subroutines shown in double lined boxes are also discussed 
individually in greater detail in the next sections. The 
subroutines shown in the hatched boxes are the adaptations 
of standard library subroutines. A brief description and 
listing of these is given in Section 10. 



Figure 8 — Block Diagram of the Compute** Progr?m 


ORIGINAL PAGE El 
OF POOR QUAIiH) 





MAIN ; In addition to the ordinary input and output files. 


the REFLCTR program (Figure 9) uses two extra files - tape 
8 and tape 20. Tape 8 is used to temporarily store a part 
of the aperture plane data while the remainder of the data is 
being sorted and rearranged. Tape 20, dictated by ^he value 
of a variable called APRDTA is used to store the aperture 
plane data both before and after being aligned along the equi- 
spaced y=constant lines. This information on Tape 20 is 
later printed out and is also used by subroutine APRPLT to 
generate plots similar to those shown in Figures 6 and 7. 

The P array contains the information on all the rays j 

emanating from the feed antenna, five quantities being 
associated with each ray viz.. Eg, E^, phase, 0', and <>’. 

MAXPTS, which is the maximum number of rays that can be 
stored in the P array, is declared in the DATA statement and is 
also used in dimensioning the P array. 

HPtlT : The purpose of this subroutine (Figure 10) is to read the 

input parameters and the feed pattern information. The 
input parameters are read from six cards which are described 
below. The feed pattern information, however, is read in a 
separate subroutine called FILL which is called by NPUT. All 
input- parameters are in the same units of measure. 

Cards 1 and 2 ; These are called title cards and contain 
any alphanumeric information in columns 1-80 for title or 
identification purposes. 
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Figure 9 *- The Main Computer Program 


SUBROUTINE NPIJT(P) 

COMMON/PARANlS/TIT'.Edh) , AORORF , XLAM,GR I D, SURFACE , APRDTA , FEED ( ' ) 
ALPHA, BETA, GAMMA, XC,YC,ZC,HFMAEX,HFMI EX, BMTP,8MPP 
NT ,NP ,NPOI NT »MAXPTS»BELLP 
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Figure 10 — Subroutine NPUT -- Continued 


Card 3 i General Information about the reflector 


AORORF 

Semi major axis of the ellipsoidal 
reflector surface, or the radius of 
the reflector sphere, or the focal 
length of the paraboloid, 

XLAM 

Wavelength in free space. 

GRID 

Spacing between the y*constant lines 
in the aperture plane. 

SURFACE 

A number that is negative, zero, or 
positive for ellipsoidal, spherical, 
or paraboloidal surface, respectively. 

APRDTA 

A variable that should be positive if 
aperture plane plots and a printout 
of the aperture plane data is desired. 

BELLP 

Semi minor axis of the ellipsoidal 
reflector surface. Even if the reflector 
is not an ellipsoid, a number should 
still be read. 


Card 4 : Feed information. 


FEED A one-dimensional array containing 

coordinates of the (x,y,z) or?.gin v;ith 
respect to the feed coordinate system, 

ALPHA, BETA, Three rotation angles that define the 

GAMMA pointing angle of the feed. (See 

Figure 2 for definition) . 

C ard 5 ; Aperture plane information. 

XC,YC,ZC x-,y-,z-coordinates of the center of 

the aperture plane ellipse, 

HFMAEX , HFMIEX Semi major and minor axes (along y- 

and z-directions) the aperture plane 
ellipse. 


YCBL, ZCBL y- and z-coordinates of the center of 

the feed shadow ellipse in the aperture 
plane. 


HFMABL , HFMIBL Semi major and minor axes (along y- 

and z-directions) of the feed shadow 
ellipse on the aperture plane. 
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Card 6 i Illumination information. 

PHI A one-dimensional array containing the 

initial, the final, and the increment 
values of $ ' in degrees over which the 
feed rays illuminate' the reflector in 
* the feed coordinate system, 

THETA A one-dimensional array containing the 

initial, the final, and the increment 
values of 6* in degrees over which the 
feed rays illuminate the reflector in 
the feed coordinate system. 

* , 

Next, the number of 6 1 values NT and the number of 

values NP are calculated. These are maximum values of 
indices i and j respectively in Figure 4. The product of 
NT and NP called NPOINT is the total number of rays emanating 
from the feed and should be less than the v?lue of MAXPTS 
which is set in the main program. 

Finally, NPUT calls subroutine FILL where E£, E^, phase, 
0', and values associated with each ray are read into the 
P matrix. 

FILL : Associated with each of the rays emanating from the 

feed antenna, there are five quantities viz.. Eg, E^, phase, 
9', and In this subroutine these quantities are stored 

in the P matrix for all the rays emanating from the feed 
antenna. One of the requirements on filling the P array is 
that the beam maximum of the feed pattern should be along the 
negative x'-axis of the feed coordinate system (Figure 1) . 

If the feed pattern is known analytically as a function of 
0 '• , and <J>' then filling the P matrix is merely a matter of 
evaluating E^(0',<J>*), E^O' ,$') , and phase for the needed 
values of 0< and <{>' and storing them in the appropriate 



locations as shown in Figure 4. For example, in the listing 
of the subroutine FILL presented in Figure 11, E£, and E^ are 
computed for each (6',$') and then stored in the matrix P 
for a diagonal horn with d«3.9 cm at 11.2 GHz. [2] . The 
phase for each ray is assumed to be zero. 

However, if the feed pattern is known in only one plane, 
or even two orthogonal planes, a scheme has to be devised to 
evaluate the E^ and the E^ components for 0£<9 *£0^ and 
g> | 1 jJjp • The subroutine FILL will have to be rewritten 
accordingly. In the next Section an example is presented 
where a circularly symmetric feed pattern is known in only 
cne plane. A simple algorithm is used to evaluate the field 
values for the other angles needed. This algorithm is shown 
in the listing of the subroutine FILL in Section 5. 

A PEPTUR : This subroutine computes the points of intersection 

of the aperture plane with the rays reflected from the 
reflector surface. If a point does not fall inside or on the 
aperture plane ellipse, its information is discarded. At 
the same time a separate count is kept of the number of points 
which lie inside and on the aperture plane ellipse. These 
are called NINTR and NEDGE respectively, and are also printed 
out. Simultaneously, the direction of the reflected ray 
corresponding to the incident ray with maximum field intensity 
is determined. In general, the direction of the beam maximum 
of the total secondary pattern is almost coincident with this 
reflected ray. The components of a unit vector along this 
ray are also printed out. If the value of APRDTA is positive, 
the data associated with the edge points and with the interior 
points are written on tape 20. 
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SUBROUTINE FILL! P »NTX »NPX ) 
COMMON/PATTERN/PHI (3) , THETA! 3) 
COMMON/MATH/PI * P 1 2 , PI 02 ,OT0R ,RTOO 
DIMENSION P(5,NTX,NPX) 

CALL SETM(0.0,P, 5*NPX*MTX) 

THR=THETA< 1 )*OTOR 

DO 300 1=1, NTX 

FAC2=PI*3.9*COS<THR)/2.6B 

FACR=COS( FAC2 ) /( 1 .0-4.0#FAC2*FAC2/ ( Pt*PI ) ) 

IF (FAC2.NE.0.0) GO TO SO 

FACC«1 .0 

GO TO 51 

50 FACC=SIN( FAC2I/FAC2 

51 CONTINUE 

PHR= ( PH I ( n-150.0)#nTnR 
on 200 J=1,NPX 

FAC1 =PI#3.9#S I N(THR)#S I M( PHR)/2.6R 
FACO=COS( FAC1 )/( 1.0-A.O*FAC1*FAC1/(PI*PI ) ) 
IF (FACl.NF.O.O) GO TO 60 
C ACA=1 .0 
GO TO 61 

60 FACA = S!NI FAC 1 )/FACl 

61 CONTINUE 

X17. = FACA#FACR 
XL v =FACC*FACO 
XMUI.1=SIN( PHR )#COS< THR ) 

XMUL?=SIN(THR)+CnS( PHR) 

Pll, I t .n*XL7.*XMllL 1-XLY*XMIJL2 
P( 2 ,1 ,.n=XI.Y«XM0Ll-XL7.*XMUI. 2 
P< 1, I , J) = AttS<P< 1 , I , J ) ) /2. 0 
P(?,I,.l)=AHS(P(2,I,J))/2.0 
P(4, I , .) ) =THR 
P<5,I ,.I)-PHR + PT 
PHR=OHR+PHI ( 3 ) *DTOR 
200 CONTINUE 

THR=THR+THETA( 3 )*OTOR 
300 CONTINUE 
RETURN 
END 
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QUANTIZ i The points of in^^section of the reflected rays 

with the aperture plane (which from now on will be referred 

% 

to simpjLy as 'points' J.n the aperture plane) are generally 

I 

not uniformly distributed. To simplify performing the double 
integration over the aperture plane, this subroutine alignes 
the points in the aperture plane along the equispaced /-constant 
lines. This aligning of points (also called quantizing of 
points) is done in two steps - first to the points on the 
perimeter of the aperture plane ellipse or to the edge points, 
and then to the internal points. The information related to 
the equispaced lines or grid bars and the resulting number 
of quantized points in printed out in appropriate format. 

P ADPA T i Performing a double integration over the aperture 
plane points, this subroutine calculates the far field 
values of E v and T. % (E being zero in the far field) for any 
0 and A. The subroutine is set up such that a far field 
pattern is requested by first picking a plane (constant 0-, 
or constant 4>-plnno) and then specifying in that plane the 
initial value, the final value, and the increment size of 
the other variable angle (<?> or 0, respectively). The output 
of this subroutine is as follows: 


A table listing the variable angle and the 

quantities E /(E). F. /(E ) . E /(F, ) , 

■* z z max y z max z y max 

* E*)/((E )* (E_)* ] in 

y y max v z y max z max 

decibels for each angle. 


2. Values of (E ) and (E ) , and 

z max y max 

3. T.ine printer plots of normalized E and E in decibels. 

z V 


A -unique feature of this Subroutine which is illustrated 
in examples 1, 2 , and 3 below is that in addition to requesting 

.far field patterns in planes and for angles specified by 

* 

numerical values (example 1) , the patterns can. also be requested 
by using 8 and $ values related to the beam maximum position of 
the secondary pattern as shown in examples 2 and 3. Refer to 
format statement number 110 in the listing of subroutine 
RADPAT in Section 8. 

Ex. 1: THETA 90.0 PHI -10.0 0.0 0.5 

III II! I | I 

I I 

I no: ooononooi o| eaaaaeiooaaoi a; aoaioni a: a o o a a a a ( aaaoooafianaennnn 

i i l a i l i ■ a ii ii mi i< i) ii !Mi ii n n ti n la » it :i n » it ti » n M n » » m n at ai a ai ai u « aa u at m ii u u m um ti u » n |i i! n aa 

In 0 = 90° plane, the far field is calculated for <f> 
between -10 to 0° in steps of 0.5°. 

Ex. 2: ,1AX-T -4.0 PHI -5.0 2.0 1.0 

I ! II 

II II I 

a a| a| no non a coo! ocooooonoooeoa aoeesoooo ooeoooooe ooooono ooonppp 1 

I I | 4 \ • 1 t I II H 111 ! M >S II HU It 19 71 77 77 74 77 71 f! 71 M M II 17 B U U » V 71 J 1 4 4 M 7 41 44 U 4 41 4 (I U SI II U M H SI SMS St I) ll I! R (4 

In 0 = ^ Q bmax” 4 °^ P lane ' the far field is calculated for 

4> values between ^bhiax~ 5 °^ to ^BMAX +2 °^ steps of 1*°°* 

Ex. 3: rtAX-P 1 0. 0MAX-T 8.0 10.0 0.1 

I II I I I 

•iiiii 

n o ; on anno on a' nj oai o| oooooocoo o o o o o o o ' o a ao a on o' on do a naan ono a a oa 

I ! I 4 S I I t I II II 17 IJ 14 IS II IMI I) II 71 77 71 74 77 71 77 71 It 3) ]l 37 r H 11 H 17 11 11 <t <1 « «} «* 41 It V « « H 11 S3 M S4 H SI 17 M SI H M 17 11 1* 

In d>= v +10°) plane, the far field is calculated for 0 

' EMAX 

values between (®bmax + 8 °^ to ^BMAX + ^°°^ ste P s of p *i°* 


APRPLT: This subroutine, a listing of which is given in 


Figure 12, is called by the main program if APRDTA is positive. 

In this subroutine the information about the points in the 
aperture plane is read from tape 20 for both before and 
after quantizing. This information was written on tape 20 
by subroutines APERTUR and QUANTIZ respectively. Two plots 
of the locations of the aperture plane points similar to 
Figures 6 and 7 are then generated but with an added feature, 
which is that different symbols are used every time the 
field strength changes in steps of 3 dB. For example, the 
output of subroutine APRPLT for a case shown in Figure 6 
will be as shown in Figure 13. 

Notice that the equivalence statement used in the main 
program makes it possible to restore the aperture plane data 

! 

from tape 20 in newly defined one-dimensional arrays j 

without requiring additional storage. 
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Figure 12 — Subroutine APR PLT — Continued 


CALL AXES( 0. »0. *90« ♦ 16. ?ZORGt SFt-1 • f 5.0*6HZ AXIS»HGT»6) 

IP=0 

I F ( ICOUNT.EQ.l ) CALL NFRAME 
I F ( ICniJNT.EQ.l ) GO TO 510 
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Figure 13 — Points of Intersection of the Reflected Rays 
and the Aperture Plane Showing the Contours 
of Constant Electric Field Amplitude 





5 . AN EXAMPLE 


A spherical reflector shown in Figure 14 is fed by a 
circular corrugated horn. The horn is 2 -polarized and its 
E-plane pattern is given. 

SOLUTION : It will be assumed that the horn has a circularly 

symmetrical far field radiation pattern with no cross polarized 
component. The FILL subroutine, which is listed in Figure 
15> is then written such that it reads the E-plane field values 
in decibels (which are also the H-plane field values) and 
calculates the field values for all the intermediate values 
of 0 ' and <|> ' . 

The geometry of the reflector is such that it subtends 
a +20*4° angle at the feed. However, as pointed out in 



Figure 14 — Dimensions of the Reflector Antenna 


Section 3, the feed pattern must also be used for angles 
greater than those needed to illuminate the reflector. In 
the present case it has been used over a 423° range. The 
boundary of the aperture plane for this axially fed reflector 
is a circle whose diameter is chosen such that it collects 
all the rays reflected by the reflector. Also# the location 
of the aperture plane is chosen such that it lies very near 

the edge of the reflector. 

* 

INPOT-OUTPUT ; A listing of input data cards and the output 
data along with the generated aperture plane plots is presented 


SUBROUTINE FILL! P,NTX,NPX> 

COMMON/ PATTERN/ PH I ( 3) »THETA ( 3) 
COMMON/MATH /PI , P I 2 , P I 02 , QTO * , R TOO 
DIMENSION P ( 5 » NTX » NPX ) 

DIMENSION E(lRl),F2(24) t IPK7) 

6320 FORMAT (5F15.10) 

CALL SETM(0.0,P,S*mJX*MTX) 

READ 6320, (E( I), 1=1, 91) 

DO 10 1=1,24 
E2( 25-1 )=F( I) 

10 CONTINUE 

DO 50 1=1,24 
11=48-1 
DO 50 J= 1 ,24 
J J=48- J 

P ( 1 , I ,J) = E2( I )+E2( J ) 

P I 1 , 1 I » J J ) =P I 1 , I , J ) 

P(1 ,I,JJ>= PI 1,I,J) 

P(1,II,.J) = P(1,I,J) 

50 CONTINUE 

DO 60 1=1, NTX 
DO 60 J=1,MPX 

PI 1 , I . J ) = 10.0** I P I 1 , I , J 1/20.0) 

PI4 t I t .J) = (THETAIl) + IT-l)*THFTA(3) )*DTnR 
P(5,I ,J) = IPHI I 1 ) + ( J-l >*PHI <3))*DTnR 
60 CONTI Ml IF 

RETURN 
END 


Figure 15 


Subroutine FILL 
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6 . SUBROUTINE APERTUR 


The purpose of subroutine APERTUR is to calculate the 
points of intersection of the reflected rays with the aper- 
ture plane and the field values at each point. 

To be able to express the components of each ray from 
the feed in the reflector coordinate system i.e., (x,y,z) 
coordinate system, the rotation matrix A is calculated. 

Then for each ray i.e., for each 0' and <*>', using the equation 
of the reflector, the following are calculated: 

R Distance along the ray between the 

feed and the reflector, 

(X0,Y0,Z0) Point on the reflector where the 

incident ray strikes, 

A one-dimensional array containing 
three components of the unit normal 
vector at (X0,Y0,Z0), 

A one-dimensional array containing 
three components of a vector along 
the reflected ray, 

A one-dimensional ar^ay containing 
three components of E along the 
reflected ray, 

Fo.int of intersection of the 
reflected ray with the aperture plane. 

Distance along the reflected ray 
between the reflector and the aperture 
plane , and 

(R+D)2tt/X + the initial phase of the 
ray . 

Y,Z,ER(2), ER (3) , and PHASE are then temporarily stored 
in array PHEW and Y,Z 0.0, 0.0, 0.0 in array PBLK . Next the 
values of Y and Z are tested to determine whether the point 
of intersection is outside, on, or inside the aperture plane 


nhat 

SR 

ER 

(Y,Z) 

D 

PHASE 
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ellipse and whether it is inside the shadow region of the feed 
on the aperture plane. Based on the results of these tests 
Y,Z,ER(2), ER(3) , and PHASE are stored in the P array with the 
internal points at the top and the edge points at the bottom. 
The same information but with a different order is also 
written on tape 20. The flow chart in Figure 16 shows the 
logic of this sorting process. Figure 17 is a simplified 
pictorial representation of the P array which shows the order 
in which the data are stored. 


If blocked set 
E * O /O and label 
it 'Blocked' 


Inside 


✓Current \ 
'point w.r.t 
^perture pi 
\ellipse > 


Edge 


1. If AFRDTA * 0, write I 
RJEW or PBLK on tape I 
20 and put appropriate! 
label* 

2. Store RiEW or PBLK 
starting at the top 
in P array. 


Outside 


If blocked set 
E ■ O/O and label 
it 'Blocked' 


1. If AIBDTA •» 0, write 
FNEW or PBLK or PINT 
on tape 20 and put 
appropriate label. 

2. Store PNEW or PBLK or 
PINT starting at the 
second position from 
bottom of P array. 


BhT S(7 
BMP >* 

BMSRX = PR(ll 
BMS3Y ■ SR(21 
BMSRZ - SR(3J 


✓ Is the X 
/present \,0\(b) X^ 
closer to Beam Max 
than the last 


Go to Next 
point 




Partial Flow Chart of Subroutine APERTUR 
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NHAT ( 3 ) =— ZO/AHRORF 
GO TO 138 

136 NHAT ( 1 )=2.0*A0R0RF/SQRT{4.0*AnRnRr**2+Y0**2+Z0**2) 
NH A T ( 2 ) = -Y0/SQRT(4.0*AnR0RF*#2+Y0**2+Z0**2) 

NHAT ( 3 ) = -ZO / SORT ( 4 .0*A0RnRF**2+Y0**2+Z0**2 ) 
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T. SUBROUTINE QUANT I Z 

The edge points and the internal points are quantised 
in this subroutine along the equispaced y«constant lines. 

First# the internal points are temporarily stored on tape 8 
and tiie edge points are quantized and then later, the internal 
points are quantized. 

EDGE POINTS : A call to subroutine FQKSORT arranges all the. 

edge points as obtained from APERTUR in ascending order first , 
with respect to the y-coord.fnate and then with respect to the 
s-coordinate and stores them in the P array starting at location 
NEDGE1 and ending at location NEDGE2 (Figure 17) . If there 
is not an edge point at A in Figure 18# a point is added at 
A and stored at the beginning of the sequence and INLO is set to 1. 



Figure 18 — Edge Points 


Similarly, if there is no edge point at B, an edge point is 
created at B and stored at the enfl of the existing sequence 

where storage locations were previously left vacant and 

* # # 

NHl'is. set to 1. NEDGEl, NEDGE2 , and NEDGE are appropriately 
modified so that NEDGEl is still the location of the first 
edge point and NEDGE2 is the location of the last edge point. | 
Next, NLO (The number of edge points on and below ) and 
NHZ (tiie number of edge points on and above Y^ ) are determined.* 
Now considering P as an one -dimensional array, the data 
associated with the NLO edge points are stored in the first 
5*NL0 locations of the p array, and the data associated with 
the NHI edge points are' stored in the next 5*NHI locations 
of the P array (DO loop ending in statement number 7500) . 

Internal points stored in these locations have already been 
saved on tape 8. However, one must check to be sure that 
(NLO+NHI)< NEDGEl so that the information is not written over j 
the memory space still occupied by the edge point data. For 
ease in understanding the rearrangement of the edge point 
data in one-dimensional form, let the first 5*(NLO+NHI) locations 
of the P array be laid out as shown in Figure 19 . The data 
are assigned to these locations as shown. The locations 
enclosed by heavy boundaries could be vacant at this point 
depending upon the values of INLO and INHI . These are assigned 
the same ER y , ER 2 , and Phase values as the point next to them 
(DO loop ending in statement number 7505) . Dictated by the 

aperture plane ellipse size and the quantizing interval, the 

ORIGINAL PAGE lb 
OE POOR QUALITY 





total number of equispaced grid lines which fall within the 
aperture plane ellipse (NBARS) is calculated. The number 
of quantized edge points (NQEDGE) and the total number of 
quantized edge points (NPOINTS) are also calculated. 

Inside the DO loop ending in statement number 8500, Z, 
ERy, ER z , and Phase aie interpolated at the end points of 
each equispaced grid line by performing a spline fit through 
the edge points along the boundary of the aperture plane 
ellipse. LI and L2 are the starting points of two temporary 
one-dimensional work arrays which are contained in the P 
array. LSTART is a number such that P(l, LSTART) , beginning 
where the interpolated values at the end points of the grid 
lines are stored, is beyond the above mentioned temporary 
work arrays. These 5 *2 *N BARS quantized values are then 
moved to the beginning of the P array. Next, the 
information about the internal points which was temporarily 
stored on tape 8 is transferred to the P array starting at 
the 5*2*NBARS+1 location. 

INTERNAL + EDGE POINTS t In the DO loop ending in statement 
number 8600, the y-coordinates of all the aperture plane 
points are aligned along the grid lines. All the points 
are then sorted and rearranged in ascending order with respect 
to the y-coordinate and then with respect to the z-coordinate . 
If the number of points along any grid line is less than 
three, that grid line is deleted. 

Finally, if APRDTA>0, the quantized data are written 
on tape 20 to be used later for aperture plots by subroutine 

APRPLT . 'uK.iNAL PAGE IS 
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8. SUBROUTINE RADPAT 


By performing a double integration over the quantized 
aperture plane points, this subroutine calculates the electric 
field in any far field direction (0,$). The double integration 
is done by first evaluating the line integrals along the y» 
constant grid lines and then performing an integration 
over these line integral values along the transverse direction 

i 

(Yg in Figure 18) . On each of the lines along which the 
above line integrals are performed, the electric field is 
known only at a finite number of discrete points. To more 
accurately evaluate the line integrals, the distance between 
every two quantized points on a line is subdivided in several 
parts and the electric field is linearly interpolated at 
these intermediate points. In the subroutine, this number 
of subdivisions has been called NPARTS . Also, to identify 
the start and the end of the y=constant lines, the y-coordinates 
of the first and the last quantized points are replaced by 
an identifying variable called SEN. Between statements number 
115 and 160, the pattern request informat n is sorted out. 

As explained in section 4, a pattern is requested by first 
choosing a plane defined by 0=constant or ^“constant (MAJOR- 
THETA or PHI) , the value of the angle defining the plane being 
called AMAJOR. Then in the plane selected above, the range 
and the increment size of the other variable angle (<fr or 0, 
respectively) is specified (MINOR=PHI or THETA, respectively) , 
AMINOR(l), AMINOR(2), and AMINOR(3) being the initial, the 


final, and the incremental values of the variable angle in 

ORIGINAL PAGh u 
OE POOR QUALITY 


degrees. 


Computation of the far field takes place between 
statements numbered 3450 and 5000. The integrals to be 
evaluated are 


B t - f f ER z coat-ei tky.ine 3 in» + kz 0 °89-»l dy<j2 


aperture 


and 
E.. = 


■// 

aperture 


(EEySin0+ER 2 cos0sin4>) 


, e j Ikysin0sin<|)+k2cos0-4] dydz 


where ER « ER , and $ are the y- and z-components and the 
y z 

Phase of the aperture electric field at a point (y,z) in 
the aperture plane. As pointed out earlier, the above 
integrals are evaluated in two steps - first by performing 
integration in the z-direction along the y=constant lines 
and computing ZI and YI which are functions of y, 


ZI. f er^os* e j(kzcos9 - t) d 2 

y=const . 
line 


and 


YI =• j (ER y sin0+ER 2 cos0sin<|>)e^ ^ zcos ® ^dz 

y=const . 
line 

and then computing E and E (called FLDZ and FLDY in the 

z y 

subroutine) as 


E 2 = FLDZ = / ZI e jkysin0sin<J> dy 


ORIGINAL PAGt L]J 
OE POOR QUALITY 



and 

8^ ■ PLOY « ^YI e 3 ky8inesln *dy. 

A flow chart explaining the steps in the above computations 
is shown in Figure 20. 

After statement number 5000, the computed values of 
FLDZ and FLDY are normalized and printed out in decibels. 

To avoid the underflow problem in computing decibels, whenever • 
the field strength is less than 10~^, -60 dB is used. 

Finally, subroutine PLOT is called to generate a line printer 
plot of the normalized E £ , E y , and the total power patterns. 

I 









/ 


c 

z 

u. 


C (A 
U. H 

z = 

« o 


I/I 


s 

u 


» * 








— a 

•*P* M« 







*t & 

a >■ 







— z 

c •> 







o « 

Z f* 







UJ » 

•i 14 







uj a 

z 







a H 

< c 







« z 

•• _J 







< tc 

— O 







H* 

o *-* 







C X 

o > 







ec a 

* 







CL w 

— c 







< z 

a -i 







► a 

a o 







a X 

X •» 







o * 

a i4 







4 X 

» •* 




a 

a 

a a 

LU LU 



• 


a 

a 

i- i- 

ec < 

O >- 




z 

z 

z X 

o z 

o c 




a 

a 

a a 

t/i a 

CO * 




♦ 

+ 

+ + 

► X 

w *4 




—a 

— » 


c * 

>- f4 



a 

a r-l 

IV 

rH (\| 

MUCL 

C O 



a 

I- — 

Wi 


ec tv -j 

-J ♦ 



£ 

z a 

a 

a a 


a > 



ec 

a c 

c 

c c 

» o a 

•— wo 



+ 

+ z 

z 

z z 

z > « 

a h- 


a 

ec 

ec < m 

*4 

*■* *■» < 

< » ► 

* *■ 


c 

c 

OI 

z 

Z Z h- 

_l O 00 G 

fsf 


z 

-> 

1— ”0 LU < 

< —* 

4 4 U 

X X 1- C 

o oo 


•—» 

4 

X 4 X II 

II X 

II II X 

*• ill- 

O H* 


z 

z 

a z i- — 

— a 

b- 

a < x of 

m ► 


< 

< 

X 4 X r* 

a X 

HNI 

ec z 4 t- 

W Csl 

ir. 


It 

IT II IT — 

W IT 

— — IT 

x of < 6 

iv* V 

•> 

a 

ec 

ii a ii a 

V M 

ec ec it 

e o 


o 

o 

a c a c 

O ec 

c o a 

a c « h t- 

_j *• 


z — 


COC2 

z c 

z z c 

z < » z c 

a 

+ 

«-• o 

< 

040m 

•- z 

— m2 

X * < ►* *• 

•- >- 

1- 

z • 

z 

4 Z 4 X 

Z l-H 

Z Z *- 

1- — H o IV 

a o 

z 

* o 

4 

Z 4 Z < 

4 Z 

4 4 £ 

z <o ui a c 

* * 

►—» 

of r-l 





»■ r-t tC Z »-* 

— IV 

c 

C a 

— * 

— * — tt A 

— * — 


a — * * a 

X 14 

a 

“5 re 

a 

P 

T 

T 

P 

a a 

T 

T 

T 

— a < a * 

a o 

z 

< *■ LT 

1 

1 1 1 1 

i i 

1 1 1 

J I z cv 

z » 

•* 

I UL rl 

X 

X X X X 

x x 

XXX 


< 

a 

o 

< 

a 

LU 

z 


t— Q. • *n 

*-i -j »- a 
H- < Z * 

V. 
t/> 

Z 
< 

ec 
< 
a 


2 Z 

c o 

Cf * 
CD Z 

3 O 
t/l O 


Ql X 
v. a 

Eg 

< _i 

z 

•s. j 

z 4 

c o 

Z *1 
Z C9 
O O 
O _J 


X >- 
KMC 
Z O _J 

' * ► u. 
ina ► 

w X 14 
CL LU O 
I- — J 

2 OU. O 


c 

II X 
LO LU 
Z -1 
uli a 
z z 
*-• c 
o o 


• 

O' 

O' 

O' 

tt 

z 

UJ 

4/i 


L/5 

H- 

ec 

< 

a. 

z 

r- • 

ii f-i 
to II 
i— *~ 

ec o f 
< < 
a a 
z a 


z z 

4 UJ 
_i oo 
X — 
"v Z 
IV H- 
ai 
00 


< < 

> » » 
Of o o 

c • o 
-) o o 

< r-l O' 
Z u. 

► ->- 
» tri — 
ILC<N 

U w 

ZHf-tL 

r- 4 0 

hCZIL 

Z < a — 

OILCIL 

Uo!ILh 


4 <<<< 4 < 4 << 

zxxxxxzzzx 

xxxxxxxxxx 

mtniriTvnminirmiPi 
• ••••••••• 

COOOOOOOOO 
UJ UJ UJ UJ UJ UJ LU LU UU UJ 


of a of of 
UICCCC 


z < < 
- £ r 


Of Of Of of of 

c c o c c 
-) “j z z z z z 

4 4 ***** ■-*•-* — 
Z £ Z Z Z Z Z 


z 

OlLU.ILILU.ILa 

O M M M 1 M k— I M H 


IL U. 


• of 

Of O 

O UJ “5 
Z X 4 

** z z 
z — < 

— I— II 

z o 

LU O UJ 

*-* o o 


Of Of Of O 
0 0 0 4- 
4- 

coca 

tt 4 tt • 

Of 1-1 IV IT | • 

O — — H LU 

I- Of Of of >• O 

O O O C X Of 

tt z z z < t- 

L*| m ih m X * 

aizia N 

c < < < ii —> ac 

II II II II 14 O X O 

ofofofofx a Q- -J 

OCr*U4XOO 
UIJlHlhO It 
QOO Oa 2 -JO 


o 

c 


in 


o 

4 


I F ( MAJOR .NE. 3 HPHI) LOOPHI-. FALSE. 

IF ( LOOPHI ) 400,3400 

COSP»COS(DEGR> PHI CUT 




a 


X 

z 


a 

UJ 

t- 

i- 1 

K 

QC 

Of ISl 

O 

< 

« 1- «• 

# 

a 

a ef sf 


»- of 

Of < isl 

a 

of # 

«• Q. ► 

to 

< — 

— of o 

K 

a. > 

isl # • 

O 

Of O' 

of — o a 

# 

* uu 

Uj I — Z; fsj 


of of of 

O o O 

tU «• •• Ul UJ 

O'O O in O Q O Q 
Z IOZ4MZMZ 

h OhIOOh Oh 

to O 00 o to o to 

H II H O II II II II 
0L h h h (L a 1“ 1“ 
Z to Z (flZWZ 

►h O •— O O i— o >-* 
tootoootooto 


If) 

« 
t- a 
to z 

O M 
o to 
* «■ 

X z 

< < H (VI 

-J -J II II 
M N C 3 

II It > UJ 

^ 5 iDZ 

ISl > H. H- 


O • 

w O 
II » 

>- • 

O Z O U: 
_i uj — r) 

LL tO II Z 
II II ►- H- 
NO>-h 

O _l II z 

JOhC 
U. >• N O 


ISl <— 

I 3 

»> >•» III 

o ~ * z 

JCUlM 
O JZ * 

M O H |(| 

•* ► (\i a. 

~ iTl 

a — a ii 
II o. — >- 
M II II Of 

o; x M uj 

UJdOC 


CTEMP*CEXP(CMPLX(0.0*ZK#Z-PH) ) 

CZ2=ERZ*C0SP*CTEMP 

CY2=(ERY*SINT+ERZ*CTSP)*CTEMP 

TSZ=TSZ+CZ1+CZ2 

TSY*TSY+CYI+CY2 



\ 


c 


c 


II II 





o 

c 





c 







o 

>- 





>- 







4 

4 





4 







4 

X 





X 








>- 





>■ 







c 

♦ 












H- 

•— 1 — H- O 





o 


>- 

> 





cc at a • 





• 


c 

C 




c 

< < < o 



t/> 


o 


4 

4 




c 

a a o. — 



~ H 




in 

m 

•Sj M 




OC X 0£ X 



o ac 


X 


• 

• 

o c 




4 4 4 _» 0. 

a. 

• < 


MwjaaNN 

^4 

w 

* 4 



~ z 

<*»> « •»*» Q_ 

£ 

z 

o a 


rsl >- a 

£ Z rsl > 

4 

4 

in m 



C 111 

c c c z 

ILr 

UJ 

► z 


c c z 

it uj u u 

rsl 

>• 

* • 



-1 to 

J J Jljl-I-O ► >- 

+ + o 

r- I- + + 

to 

to 

M w 

H H 


c • 

OOC~* 

o o 

• ^4 C 

c c — 

O O «H M 

K 

r- 

4 4 

+ + 

o 

— o 

Q. 

4 

4 

O II 

+ 

j j a 

4 4 isl >• 

+ 

+ 

•si V 

osin 

» UJ 

rsl >- I X 

C 

o 

— z o 

C O x 

O O o o 

rsl 

>- 3 

UJ to l/) 

_i uj 4 

UJ H t 

1 1 3a 

-1 

_l 

II 

— 1 

M M UJ 

-JJ+ + 

UJ O 

c UJ UJ 

3 h- i- 

o z 


3 — O 

HHlLJ 

w 

O > O O N >• O 

OON>NN 

3 _J 

-J 3 Z 

z + + 

M M 


2a J 

rsl >- Z II 

►— 

*— < 

to O > 

II II II 

MMtOt/SfSl>- 

Z U. 

UL Z >- 

* »— ►— « 

II II 

o 

II o 

— — • a 

rsl 

>- *- m 

II 

O Q a 

NH-I-OO 

M II 

II M || 

t- ISJ V 

O 3 1- H 2 >■ 

II II — z 

II 

ii 

II 4 

O 

J J£ 

II II II II II 1 

i— M 

> h- C 

Z II II 

-1 UJ 


Z UJ — 

MM It UJ 



fsl 

_J 

C O UJ 

NNmVhhZO 

C z -J 

o ~ <~ 

o z 

CCZU. 

M >- > y- 

N >■ 1/1 C O 

M M J- 

M > tO tO M >- 

D —1 

-ICC 

O ISI V 

M — « 

o o >• •— * 

ooaoooi— o>- 

fsJ>-UOOI-l-OOOU_ 

UL l_) >- 

O 

o 


o 

o 







o 

o 

o 

o 


o 

o 







O 

o 

h- 

O'. 


o 

(M 







m 

4 

m 

fO 


4 

4 







4 

4 



V> 

c 

2 

< 


e 

z 

uu 


o 

o 

O' 


rr 

PM 





M~ 




c 

I 

X 



K 

1 - 

H 




z 

Z 

o 


P 

mm 

w 

o 


O 

Ml 

>- 

o 



o 

c 

in 


p» p> - — 

_J 

-J 



Z Ml >■ 

X 

X 

o 


X- c c 

MM 

MM 

H- O 

l/l _J —1 

X 

X 


c 

X x 

M 

M 

p 

p* 

• mm ppp 

Ml 

>- 

o 

m 

X If) to 

X 

X 


► 

— z c « 

«3 

< 

pm 

o 

o •< < 

Z 

z 

Of 

o 

• o u x 

X 

M 4 

X 

O — II It 


MM 

X 


*• 3 .-. MM 


pH 

p 

PM 

O X I X 

X 

X 


MM 

• Z 1 - H < 

< 

* 

X 

o *-« z z 

z 

z 

Of 1 - 

a 

— ♦ w W 

< 

< o o 

P 


« 

CD 

c 

DC 




2 


MPP 





X 


o 





ft 


MH 



MM 


#• 


e 



o 


X 


c 



M 


x 





o 




< 



o 


« 


ft 



_» 


>- 


c 



<. 


V. 


• 



ft 


>■ 


— c 



o 


X 


^ mM 



• 


c 


M- II 



o 


ft 


«P -M 



M* 


► 

Ml > M. 



II 


X 

c 

C — 



— — Cd 


X 

-J 

-1 Of 



Ml >- O 


•> 

X' 

X X 



X X z 



MM 

M- 3 



< « Of 


> 

X 

x a 



Z £ 3 


\ 

PPM 

MM 



x x a 


Ml 

O 

o 



w W 


X 

pH 

mM c 



o o — 


c 

o 

c -> 


X 

P-4 M c 


ft 

c 

C 1 


o 

C C <-■ of 


*- 

_• 

_J X' 


of 

PCI c- 


X 

< 

< a 


2 

J U' 2 


>0 

ft 

ft • 


a 

< < o •- 

(f'. 

p 

c 

C -M 



ft * • z 

• 

ft 

• 

• • 


> 

O O M < 

ec 

Ml 

o 

C H- 


>- 

• • • * 

X 

M. 

IM 

X O 


X 

O o 1 - of 

fO 

>- 

II 

II • 

PM 

o 

l\! fM C C 

p 

£ 

Ml 

V IV 

K 

•> 

ii n • — z 

in 

C 

fiC 

X ft 

ft 

> 

scco>->-m- 

< 

ft 

c 

C ft 

ft 

Ml 

o o X X z 

p 

*• 


MM 

MM 

X 

Ml >• « < ► 

X 

X 


PP mm 

PH 

c 

Z Z Z Z of 

in 

in 

o 

O MPP 

MM 

*■ 

X X X X c 

-> 



— > 

V 

Ml 

ft ft "3 

m 

ft 

C 1 

1 c 

o 

>- 


t 

Ml 

• X' 

X _> _l 

X X 

c o x x z 

cc 

N 

o o 

c a. 

IU 

c o 

H M < < <1 

X 

Ml 

pO • 

t *«hi 


T * 


|| || II <-4 Ml >• II II m* • C 

OOh-DQN>-0 _1 

_i_iva-j-jxx + c- 

OO II — XI X « < O — 

•-•MMU.t-ii-'XZ ii u. a. 

M>-M>->lLULU.U.O^*- 


X of 

o a 

z z 

P-M P-M 

V- T 

z < 

O II 

o o 


•o 

I 

II 

u 

C I I 

£ll‘U! 

of o o 

3 • • 

a. -i ■-! x x < 
ii • • < < z 
SC *— I— 

coo 

X Ml > 

a. x x 

ii < < 

£ Z Z 

C X X 

Ml — 

Z X X 

IL *-M I— I 


x x ► ► as 

u. x re in & c 

+ o < c # 

Ml Ml -3 - Z - X 

X M« O 

H J pJ 

z x •• < 

U_ IL » <-m * w 

ft ft O I X X 

N N O >-l «0 

x x x — x — 

<3 •< I— I— 

Z Z K- 4 t- « 

x x z z z z 

— 

C of c 
x a x 


x Of 

*-i Q. 


I <-4 u. 

M • • + 

— I- H- (V 
— O O * 
I- — • • ft 

Z if ~ 

•> yj m *— * ** 

mh 3 — ~ — 

ii a m 

M II O 

v _) 

O S 3 UJ 

o c »- 

r- ii a 

Ml — 

2 t 11 . 

o c •— 


o c 

.J -J 

a & 

or a 

c o 

IL IL 


>- Ml 
C O 
— • _J 
X Up’ 

M N« 

u. x 

•p* 

u. a. 


x . . of mi • 

+ X MJ t-i 

(v a eo <-* • 

♦ « SIS I ox 

ft C c C C »-in-»Nl>- 

* N > >■ Mm O » Ift X X 

— X X X x «n — C c 

— X X X X of » •OCII II 

Ml I I I I X O O' o — — > 

OM>N>XO'XZmm 

JIIICHO-hp-- 

XCCCOII X Z Ml > 
M II || II tt CD X « < O O 
XIM>->-M czx + -> _l 

M>N>tfM(fOUJUJ 

cc x £ i 3 a; C ii mm 
ococaaxoxx 


ss* 


-SSS 


o 

c 

o 

pO 

o 

c 

c 

c 

O' 

* in 

<0 

<c 

X 




Ift 




— 

*■ 

» 


ft 

• 






® 

s 



• 

o 

o 



• 

*■* 

w 



• 

• 

z 

z 


« 

• 

oc 

oc 



• 

Ul 

UJ 


» — * 

• 

H 

I- 


r- r- 

in 

N- 

r- 

— 

• • 

*4 

4 

« 

ffi 

in 


a. 

a 

o 

—4 H 

z 



— 

o & 

c 

>- 

>- 


« * 

•4 

of 

a 

z 

ft ft 

N- 

4 

4 

Of 

N H 

4 

C 

C 

u. 


e £ 

Z 

z 

N- 

ft ft 

c 

c 

c 

N- 

* 

IU 

o 

o 

4 

N- N- 

N> 

Uj 

u. 

a 

• • 

Z 

10 

10 


lf> m 

•4 



oc 

— ( H 


u. 

u. 

w 

IS u 

Of 

c 

c 

3 

•> - 

C 



C 

ft ft 

u. 

N- 

K 

a 



Z 

z 


ffl c o 

c 

UJ 

IL' 

c 

CC J 

UJ 

Z 

z 

U' 

> _» -J 

IO 

c- 

c 

ISl 

ZOO 

D 

CL 

& 

— 

U.NM 


z 

X 

_> 

•> N R 

of 

c 

c 

4 

> ~ -* 

IM 

o 

o 

Z 

X — ~ 

& 

1 

1 

OC 

4 M >- 

z 

— M 

>- 

c 

Z 1 1 

3 

CM 


z 

u. o e 

z 

• c 

c 


^ j j 


<C UL‘ 

UJ 


eoujiu 

z 

U. M 

NJ 


O M M 

o 

* 14 

>4 


NU.U. 

— o 

O -1 

C -J 

c 

Z — — 

NUJ<< 

U' 4 

u. 

u. X X 

< c 

• ? 

c z 

c 

* 4 4 «/> -1 * 

X O' — 

»af- 



N * * *- O or O o I of C X of — 
X-»d£ftCNZ(-CZhC3 X 
"5 * Z ”0 Z "0 K 

zoCaiMvi »<x »<xz 

u. J JZI-£SONIO>-lO * 

o o, z vine in o ino( 

»N IM • •— -ft I ► — _J • w Ui 

o it ifi-^»-»L»nh-ujirK*o 

uj in * ft in « n ttC^-r-C<-i^CQ.co.' 
Dr-— r-N-ju.Nju.hj ► -* d 

z n- n- ho.* a* a® zz 

— t- < I- 4 I- < « N- EC N- O C — Of 

N-ZZ ZZZZ_tDZ_l2Z-Jrr-»-D 

Z t-OC-K JN« JNw zr-o 

C of C C< JOOIUZ 

o a u. bu.uiuaaou.aoaoou!uu 

• • • • 


t 


c 

o 

in 

in 

c 

in 

in 

N- 

r- 

r- 

r- 

N- 


o 

c 

c 

o> 


9. SUBROUTINE INTERP 


INTERP is called by APERTUR to interpolate a point at 
the edge o£ the aperture plane ellipse when, in the sequence 
of sorting out the reflected rays, the point of intersection 
of a reflected ray with the aperture plane falling inside 
(outside) the aperture plane ellipse is followed by another 
which falls outside (inside) the aperture plane ellipse, the 
two points not being on the opposite sides of as shown in . 
Figure 21. 

Of these two points, based upon ENEARl and ENEAR2, 
the point closer to the ellipse circumference is selected. 

The y-coordinate of this closer point, YEDGE is then taken 
as the y-coordinate of the interpolated edge point and a 



Point, #2 


Figure 21 — Interpolated Point at the Edge of Aperture Plane Ellipse 


corresponding value of z-coordinate , ZEDGE is calculated 

using the equation of the ellipse. Next, calculating the 

distances Dl and D2, £R . ER , and Phase are linearly 

y * 

interpolated at the edge point and assigned to PINT (3), 
PINT (4) , and PINT (5) , respectively. 





SUBROUTINE I NTERP( PI »P2 ,PINT ) 

COMMON/PARAMS/TITLEI 16) t AORORFtXLAMtGR 10* SURFACE tAPRDTA* FEED! 3) * 
ALPHA, BETA » GAMMA tXC»YCt2C,HFMAEX,HFMIEX*8MTP*BMPP* 
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CONTINUE 

RETURN 

END 


10. OTHER SUBROUTINES 


MULT32 : This subroutine multiplies a 3x3 matrix ORA by a 

3x2 matrix OPB to yield a 3x2 matrix ROP. A listing of the 
subroutine is shown in Figure 22. 


SUBROUTINE MULT32( ROP ,OPA,OPB ) 

DIMENSION R0P(3,2),0PA(3,3),0PB(3,2) 

ROP (1,2) =OPA( I » I )#OPB( 1*2 )+OPA ( 1 ,2 )*0PB( 2,2)+OPA 11*3 )*OPB( 3,2 I 
ROP ( 2 • 2 ) *=OPA( 2.1) #OP B { 1 , 2 ) +0P A ( 2,2) *OPB (2,2) +OPA (2,3) *OPB (3,2) 
ROP ( 3 , 2 ) »OPA ( 3 , 1 ) *OPB ( 1 , 2 ) +OP A ( 3 ♦ 2 ) «nPB ( 2 , 2 ) +OPA ( 3 ,3 ) *OPB ( 3 , 2 ) 
R0P(1,1 )=0PA( 1,1 ) ^OPS ( 1, 1)+0PA( 1,2)*0PB( 2,1 )+0PA( l,3)*OPB(3, 1 ) 
ROP (2,1 )=OPA(2,l )*OPB( l,l)+OPA(2,2)*OPB( 2 ♦ 1 )+OPA { 2 ,3 )*OPB( 3 , 1 ) 
ROP( 3,1 )=OPA{ 3, 1 )#OPB( 1,1 )+OPA ( 3,2)#OPB(2,l )+OPA ( 3, 3 )*OPB( 3, 1 ) 
RETURN 
END 


Figure 22 — Subroutine MULT32 

FQKSORT : This subroutine arranges in ascending order the NN 

sets of numbers in the A(LTH,NN) array according to numbers 
stored with LTH=1,2, . . .etc. A listing of this subroutine is 
presented in Figure 23. 

CURVl t This subroutine determines the parameters necessary 
to compute an interpolating spline under tension through a 
sequence of functional values. The slopes at the two ends 
of the curve may be specified or omitted. For actual compu- 
tation of points on the curve, it is necessary to call the 
function CURV2. At the time of input: 

N Number of values to be interpolated, 

X An array of the N increasing abscissae 

of the functional values, 
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Figure 23 — Subroutine FQKSORT 
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igure 23 - Subroutine FQKSORT — Continued 


I • 

Y An array of the N ordinates, i.e., Y(K) 

is the functional value corresponding to 
X(K) , 

SLP1 and The desired values for the first 
SLPN derivative of the curve at X(l) and 
X(N), respectively. If the quantity 
SIGMA is negative, these values will be 
determined internally and the user need 
only furnish place holding parameters 
for SLP1 and SLP2 , 

YP An array of length N, 

TEMP An array of length N used for scratch, 

SIGMA Contains the tension factor. This is 

non zero and indicates the desired 
curviness of the spline fit. If the 
absolute value of SIGMA is nearly zero, 
(e.g. 0.001), the resulting curve is 
approximately a cubic spline. If the 
absolute value of SIGMA is large (e.g. 
50.0), the resulting curve is nearly 
a polygonal line. 

| Upon return from this subroutine YP contains values 

i 

§- 

I proportional to the second derivative of the curve at the 

| given nodes. N, X, Y, SLP1, SLPN, and SIGMA are unaltered. 

A listing of this subroutine is presented in Figure 24. 

CURV2: This function subroutine interpolates a curve at a 

given point using a spline under tension. The subroutine CURVl 
should be called earlier to determine certain necessary 
parameters. At the time of input T contains a real value 
to be' mapped onto the interpolating curve , N , X , Y , YP , and 
SIGMA are defined the same as in CURVl, and IT is an integer 
switch. If IT is not 1, this indicates that the function 
has been called previously with the same N, X, Y, YP, and 
SIGMA but the current value of T exceeds the previous value. 
Upon return, CURV2 contains the interpolated value. A 
listing of this function subroutine is shown in Figure 25. 




SUBROUTINE CURV1 ( N* X , Y, SLP1 . SLPN, YP» TEMP,S IGMA ) 
INTEGER N 

REAL X(N) ,Y(N), SLP 1 ♦ SLPN, YP( N ) ♦TEMPI N), SIGMA 
NM1 * N-l 
NP1 = N+l 
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Figure 24 - Subroutine CURV1 


DX1 = DX2 
UIAG1 = 0IAG2 

DIAGIN = 1 . / ( DI AG1—SPDI AG* TEMP ( NMl ) ) 

YP ( N ) = DIAGIN*(SLPPN-DX2-SP0IAG*YP(NM1) ) 
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Fiours 24 — Subroutine CURV1 — Continued 


FUNCTION CURV2 ( T.N. X., Y. YP. SI GMA , T T ) 

INTEGER N.I T 

REAL T * X ( N ) * V I N ) » YP ( N ) *SIGMA 
S = X(N)-X( 1 ) 

SIGMAP » ABS(SIGMA)«FLOAT(N-l)/S 
IF UT.EQ.l) 1 1 ■ 2 

1 DO 2 I = 1 1 » N 

IF IXIII-T) 2.2.3 

2 CONTINUE 

I * N 

3 IF (XU-1) .LE. T .OR. T .LE. X(l)) GO TO 4 

II = 2 
GO TO 1 

4 DELI = T-X(I-l) 

DEL2 = X ( I ) -T 
DELS = XU)-X(I-l) 

EXPSi = EXP ( SIGMAP40EL 1 ) 

SINHD1 * . 5*(EXPSl-l. /EXPSI ) 

EXPS = EXP(SIGMAP*OEL2) 

SINHD2 = . 54 (EXPS- 1. /EXPS) 

EXPS = EXPS14EXPS 

SINHS * • 54 ( EXPS- 1‘. /EXPS ) 

CURV2 = (YP(I )4SINHOl+YP( 1-1 )4SIMH02)/SINHS+( (Y( I )-YP( 
1 I ) )4DELl+( Y( 1-1 )-YP ( I -1 ) ) 40EL2 ) /DFLS 

II = I 
RETURN 
END 


Figure 25 — Function CURV2 


PLOT : This subroutine generates a line printer plot using 

NT values stored in the F array and labels the plot as NAME. 

FMAX is the maximum value of F. The points less than -60 dB 
are indicated by ' <* sign on the horizontal axis. The listing 
of this subroutine is presented in Figure 26. 

SETM and MOVEM: These two small subroutines are used throughout 


the program. The subroutine SETM sets the first N consecutive 


SUBROUTINE PLOT ( NAMI= * FMAX « F*NT ) 

DIMENSION NUM ( 1 2 ) * I FRM ( 4 ) ,F(1) .NAME! 51 
OATA I FRM/O » Ot 0» 6HX * *« * ) / 

DATA NUM/IHO » 1H1 » IH2 » 1H3* 1H4» IH5 * 1H6* 1H7 * 1H8 » 1H9 »2H10»2H1 1 / 
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Figure 25 ~ Subroutine PLOT 


locations of the array ANEW to x. The subroutine MOVEM is a 
little different. It transfers the first N consecutive 
numbers from the array AOLD to the N consecutive locations 
in the array ANEW. A listing of these two is presented in 
Figure 27. 




SUBROUTINE MOVFMl AOLD* ANFWtN ) 
DIMENSION AOLD(l),ANEW(l) 
IFtN.LE.OJ RETURN 
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11. CONCLUSIONS AND REMARKS 


The computer program presented here has been successfully 
used to compute the radiation properties of large spherical 
reflectors [3] and to study the optimization of parameters 
for wide band radiometric applications [4j. The computation 
time is very reasonable, e.g., on a CDC 6600 computer, the 
computation time per far field point has been found to be 
0*63 sec for a 700 X diameter spherical reflector. Since 
the aperture integration time required for computation of 
the far field depends upon the number of aperture data points, 
which in turn is determined by the angular increment between 
the rays used in the feed pattern, the computation time is 
somewhat insensitive to the absolute size of the reflector. 

There are only three places in the entire program where 
a decision is made based upon the reflector surface i.e., 
whether it is paraboloidal, spherical, or ellipsoidal - once 
in the subroutine NPUT and twice in the subroutine APERTUR. 
These places have been labeled with ******** in columns 73-80. 
In the subroutine NPUT the surface information is used to 
output the information in appropriate format and in the 
subroutine APERTUR it is first used to compute R and then to 
write the three components of the unit normal to the reflector 
surface, NHAT(l), NHAT ( 2 ) and NHAT(3) . Modifying the NPUT 
and APERTUR subroutines at these three places, the program 
can be used for any other reflector surface that can be 
expressed analytically. 
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